Measuring spectrum of spin wave using vortex dynamics 
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We propose to measure the spectrum of magnetic excitation in magnetic materials using motion of vortex 
lattice driven by both ac and dc current in superconductors. When the motion of vortex lattice is resonant with 
oscillation of magnetic moments, the voltage decreases at a given current. From transport measurement, one 
can obtain frequency of the magnetic excitation with the wave number determined by vortex lattice constant. 
By changing the lattice constant through applied magnetic fields, one can obtains the spectrum of the magnetic 
excitation up to a wave vector of order 10 mrr 1 . 

PACS numbers: 74.25.Uv, 74.25.F-, 74.25.Ha 



I. INTRODUCTION 

Measurement of the spectrum of magnetic excitation is cru- 
cial for understanding magnetic properties of magnetic mate- 
rials. Neutron scattering has been widely used to measure the 
spectrum. However, the source of neutron is not easy to ac- 
cess. It also requires large samples. For example, the neutron 
scattering is not applicable to ultra-thin films because the in- 
teraction between the spin wave and neutrons is week 1 . On the 
other hand, conventional techniques such as the ferromagnetic 
resonance and Brillouin light scattering can only measure the 
energy gap of spin wave. 

Abrikosov vortices being magnetic excitation in supercon- 
ductors are expect to interact strongly with the magnetic mo- 
ments, which points a possible way to measure the magnetic 
excitation using vortex dynamics. It was proposed that one 
can use vortex dynamics driven by a dc current to measure the 
spin-wave excitation through transport measurement of the IV 
characteristics 2 . One may be able to obtain the spectrum of 
spin wave up to the wave number 1 with £ being the coher- 
ence length. The proposed technique can be applied to super- 
conductors with coexistence of magnetic and superconduct- 
ing ordei^. One can also measure the spectrum of magnetic 
excitation in conventional magnetic materials by fabricating 
artificial bilayer systems, consisting of the magnetic material 
to be measured and a superconducting layeiR 

In this work, we propose to measure the spectrum of mag- 
netic excitation in magnetic materials using motion of vortex 
lattice driven by both ac and dc current through transport mea- 
surement. The advantages with ac current are: 1) it does not 
require large current to reach the resonances between the vor- 
tex motion and spin-wave excitations, thus the nonequilibrium 
effect that leads to the instability of vortex lattice and heating 
effect can be minimized; 2) the effect of random pinning cen- 
ters can be minimized with an ac current; 3) additional infor- 
mation such as wave vector of the vortex lattice can also be 
obtained. 



II. PROPOSAL 

We consider a magnetic superconductor where magnetic or- 
dering coexists with superconductivity as shown in Fig. [TJa) 
or a bilayer system with magnetic film and superconducting 



film as shown in Fig. [TJb). We assume that the magnetic sys- 
tem has an in-plane easy axis (x axis). The spectrum of spin 
wave Q(k) is gapped due to the anisotropy, where k is the 
wave vector. A magnetic field to create vortex lattice is ap- 
plied perpendicular to the easy axis of the magnetic moments 
(z-axis). A transverse current both with dc, 1^, and ac com- 
ponent, I ac sin(w/f), then is applied to drive the vortex lattice. 
Driven by the Lorentz force, vortex lattice moves with a dc 
velocity v^ e , and also oscillates due to the ac current and the 
interaction with the magnetic moments, v = Vrf c +v flf . The mo- 
tion of vortex lattice perturbs the magnetic moments and ex- 
cites magnons. The resonance between the motion of vortex 
lattice and precession of magnetic moments is achieved un- 
der appropriate condition [see Eq. ( 14 1], which can be probed 
from IV curve. 

Using the vortex-as-particles approximation, the motion of 
vortex is described by an over-damped Langevin equation 



r)d,ri = F p + F L + F vv - d n Hnt(r - r,), 



(1) 



where r, = (xuyi) is the vortex coordinate and r\ = 
B z H C 2cr„/c 2 is the Bardeen-Stephen damping constant with 
H C 2 the upper critical field and <j„ the normal state conduc- 
tivity just above the critical temperature. F vv is the repul- 
sive force between vortices. Fl = Fd c + F ac sin(w/f) is the 
Lorentz force due to the ac and dc current. F p is the pinning 
force, and < Ki nt is the Zeeman interaction between magnetic 
moments and vortices 



-r,-) = - J" B v ^(r- 



Ti)M z (r)dr 



(2) 
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FIG. 1. (color online) (a) and (b): Schematic view of vortex dynam- 
ics in (a) magnetic superconductors (MSC) where superconducting 
and magnetic ordering coexist, (b) artificial bilayer systems consist- 
ing of superconducting (SC) and magnetic (M) layers. 
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where B vz is the magnetic field associated with a vortex and 
M, is the magnetic moment along the z-axis. We have used 
a continuum description of the magnetic subsystem because 
the vortex size is much larger than the lattice constant of spin 
subsystem. The random pinning centers distort lattice order 
in static 6 . However in the flux flow region, the vortex lattice 
order is recovered because the motion of lattice quickly aver- 
ages out the effect of random pinning centers-^! In this region, 
it is safe to approximate straight vortex line along the z-axis 
and the problem becomes two dimensional. 

In the following calculations, we will focus on the reso- 
nance between motion of vortex lattice and magnetic subsys- 
tem, and neglect the resonance due to pinning centers. We use 
an approximation that the motion of vortex lattice is not af- 
fected by the magnetic moments. The motion of vortex lattice 
in the presence of the dc and ac current thus is described by 



r,(f) = R, - y dc t - r ac sin (ajjt + dj) 



(3) 



where R, forms a regular lattice with a wave vector G and <p 
is an arbitrary phase. The vortex lattice oscillates with fre- 
quency 0)1 and amplitude r ac = J a c^()/(cria>i) due to the ac 
input current J ac . 

We use the quasistatic approximation that the structure of 
vortex driven by the Lorentz force remains the same as the 
static one. The magnetic induction B of the moving vortex 
lattice is described by the London equation taking the mag- 
netic moments M into account^^E-D 

X]y x V x (B - 4ttM) + B = O Y 6 [r - r,(f)] % (4) 

i 

with Ai the London penetration depth without magnetic sub- 
system, 3>q = he I (2e) the quantum flux and z the unit vector 
along the z-axis. 

The magnetic field "seen" by the magnetic subsystem is 
H z = B z - 4nM z . Using the linear response approximation, 
M,(k, to) = Xxfe' &0.ff z (k, co) with^fk, to) the magnetic sus- 
ceptibility, we obtain the induced magnetization M-(k, at) = 
Xzz(k, u>)B z (k, cl>)/[1 + AnXzziK to)]. We then obtain the mag- 
netic fields associated with the vortex lattice from Eq. (H| 



B z (G,to) = B(fK(G,a)) 



y\ +4n X zz(G, cS) 



+ 1 



(5) 



with Bo the average magnetic induction. The London pen- 
etration d ept h is renormalized in the presence of magnetic 
subsysterrPEEI] The function 7C(G, to) is 



K(G 



,oj)= f 
Jo 



dt exp [/G • (yd c t + r ac sin (ajjt + <p)) - icot] 



(6) 

Using the Fourier expansion of the Bessel function, we have 

n=+oo 

7C(G, to) = 2^ Jn(G- r ac ) exp(incf>)5 (G-\ dc -to + ncoi) 

n=—oo 

(7) 

where Six) is the Dirac delta function and J„ is the Bessel 
function of the first kind with an integer n. Due to the pinning 



centers, the Bragg peaks of the vortex lattice are smeared out. 
The velocity of vortex is also fluctuating. These two effects 
can be taken into account by replacing the 5 function in Eq. 
(7}by 



W{G,d) 



1 



G • Vrf e — to + nui + iy ' 



(8) 



where y accounts for the broadening of the resonance peaks 
by the pinning centers. 

The IV characteristics of the system can be derived from the 
power balance equation. The energy input per unit volume is 
JE with J the external current and E the electric field. The 
power dissipated per unit volume by qu as ip articles due to the 
motion of vortex lattice is rjv 2 . The velocity of vortex lattice 
can be measured from voltage E according to v = Ec/B^. 
The energy per unit volume P{E) transferred from the vortex 
lattice into the magnetic subsystem is^^l 



P(E) = - (M z (r, t)d t H z (r, f)> 

= - f lm\x zz (k, cj)]to \H z (-k, -co)\ 2 cfkdco 



(9) 



where (• • • ) represents time and space average. The trans- 
ferred energy P > finally dissipated through magnetic 
damping. Using H z = B z /(l + 4nx zz ), we have 



G,« 



dtolm \xzz (G, w)] oj 



BqJ„ (G • raj W(G,to) 



A 2 L G 2 + l+4n Xzz (G,to) 



(10) 

From the power balance condition JE = r/v 2 + P, we derive 
the IV characteristics 



J-^E+'-PiE). 



(11) 



In the presence of magnetic subsystem, the effective viscosity 
of vortex is enhanced 



left =T} + 



P(E) 



(12) 
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Voltage V 

FIG. 2. (color online) (a): Schematic view of the IV curve and (b) 
dl/dV curve. 
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For a given current, voltage drops because of the energy ex- 
change between the vortex lattice and magnetic subsystem. 
The susceptibility of the magnetic subsystem isP 



*zzO, k) 



Q 2 (k)-w 2 + iw/3' 



(13) 



where co M = fj. 2 n M /(2ft) with fj. the magnetic moment and n M 
the density of magnetic moment. £2(k) is the dispersion of 
spin-wave excitation and B the relaxation rate of spin wave. 
Therefore the resonance takes place when 



rw, + G(B) • \ dc = Q(G), 



(14) 



is satisfied. The principal axis of the moving vortex lattice is 
along the driving direction, and G(B) ■ v<f c = Invdc V^o/^o 
for a square lattice. The resonance amplitude depends on 
the amplitude of the ac current according to [J„ (G • r ac )] 2 . 
Both the random pinning centers and magnetic damping con- 
tribute to the broadening of resonance, yielding a linewidth 
of resonance of order B' = B + y. Away from the resonance 
|Q 2 (G) - (cl>'„) 2 \ » co' n B' with co' n = G ■ \d c + not/, we have 
Im[^ zz ] m 0. The current is then given by J — rjc 2 E/B 2 , be- 
cause there is no spin wave excitation with frequency co' n at the 
wave vector of vortex lattice G. At resonance, the current is 
enhanced for a given E, J = t]c 2 E/B 2 ) + AJ. Here we estimate 
the current enhancement AJ. 

We consider frequency slightly away from the resonant fre- 
quency Aco = O - (ncjj + G ■ \ dc ) with B' « Aco « Q. We 
then have Re [xzz] = com/ '(2Aa>) and Im [xzz] = com/3/(2Aco) 2 . 
For a square vortex lattice, we have G = 2n ^fBo]<5^) (l x , lyj 
with integers l x and l y . In the interval of frequency deviation 
Aco where A 2 L G 2 » 4ttx zz (G, co' n ), we estimate the enhance- 
ment of current over the linear background when the current 
is injected along the x direction 



AJ 

T 



16n 2 Tin M B'B l 2 x J 2 (G • r ac ) 



to' n V®o 



(15) 



For the HoNi2B2C magnetic superconducto^^EIS we have 
H c2 = 10 T, cr„ ~ lO^Q-mU 1 , n M = 10 22 cm" 3 and 
B ~ 10 6 s~'. As y depends on the concentration of pin- 
ning centers and is not known for magnetic superconductors, 
here we neglect y and only keep B for the estimation of AJ. 

We then estimate AJ/J * 0.8Z 2 7 2 (G • r ac )(/ 2 + iff 1 at a fre- 
quency co' n ~ 10 GHz. 

To achieve measurable enhancement of current AJ at reso- 
nances, one requires G • r ac ~ 1. Using r ac = J ac <f>{)/(cricoi), 
we estimate the amplitude of the ac current that yields mea- 
surable resonances 



TjiOjC _ to i B c2 cr n 
GOn ~ G c 



(16) 



For HoNi2B2C we estimate J dc ~ 10 9 A/m 2 when w/ m 
10 GHz, G ~ 0.3nm _1 , which is much smaller than the de- 
pairing current J Ap = cB c2 ^ j{6 V3^/l 2 ) ~ 10 13 A/m 2 with 
^ w 50 nm and A L m 100 nm.^Thus the measurable current 
enhancement AJ can be realized experimentally. 



Here we present a procedure to extract Q(G) from the IV 
curve. One measures IV curve at a given vortex density Bo 
and ac current. When the resonance condition Eq. 



satisfied, the current is enhanced according to Eq. ( 15 



14} is 

i for a 

given voltage, see Fig. [2^a). This enhancement can be seen 
clearly from the curve of dl/dV as a function of V, see Fig. 
|2jb). The linewidth of the resonance is due to the magnetic 
damping and random pinning centers. The voltage where cur- 
rent is enhanced is the resonant voltage. One then changes 
the frequency of the ac current coj, and measures the resonant 
voltage as a function of coj. From the voltage, one knows the 
velocity of vortices lattice at resonance. One then plots coj 
as a function of the resonant velocity. From the interception 
with the vertical axis, one obtains Q(G) according to Eq. (14) . 
From the slope, one obtains G along the driving direction. By 
changing G through external magnetic fields, one then obtains 
the spectrum of spin wave f2(G). The spectrum can also be 
determined from the sub-resonance with n-2. 



in. DISCUSSION 

Here we discuss the region that G can be tuned. The mag- 
netic field associated with vortices cants the magnetic mo- 
ments and induces magnetization along the vortex axis M z 
with M z = XzzBq/(4-ttXzz + 1)> th us the spectrum measured 
is Q.(G,M Z ). To get the spectrum at M z = 0, we need 
M z /M s <K 1 with M s the saturation magnetization, which 
gives an upper bound for Bo- Since Bo ~ ^o/a 2 with the 
vortex lattice constant a, this limits the maximal wave vector 
that can be achieved G ~ 2n/a. For the HoNi2B2C magnetic 
superconductojEEEl, M s * 1000 G and X zz ~ 0.03, which 
gives maximal G max ~ 0.3 nirT 1 . For a large susceptibility 
Xzz ~ 1, it was shown that vortices may form clusters due to 
the attraction between vortices 12 . To create vortex lattice, the 
minimal vortex density is n v » l/A 2 L . Thus the lower bound of 
G in this case is G = 2n/Ai. For^- « 1,G can be tuned to 
corresponding to a dilute vortex lattice. 

The magnetic correlation length % m is much smaller than the 
maximal vortex lattice constant, £ m G max <K 1. Therefore we 
can expand Q. as £2(G) w co g + h 2 G 2 /(2m s ) for a ferromagnet, 
with cjg the energy gap and m s the mass of spin wave. For an 
antiferromagnet, £2(G) ~ co g + \ s -G with v s the spin wave ve- 
locity. Although the present method can only measure small 
portion of the Brillouin zone compared with the neutron scat- 
tering, addition information such as the mass m s or velocity v s 
of the spin wave can be extracted, which is advantageous over 
the ferromagnetic resonance measurement. 

For magnetic materials, typically the energy gap is co g ~ 
10 GHz to 100 GHz. Experimentally, one can achieve coj ~ 
100 GHz by microwave. Thus one can measure the spin- wave 
spectrum at a low velocity of the vortex lattice. The nonequi- 
librium effect that destabilizes the vortex lattice, caused by the 
Larkin-Ovchinnikov mechanism 16 , can be avoided. 

The ac current itself generates ac magnetic field Z/ ac = 
2nJ. dc Llc with L being the linear size of the system. The in- 
duced magnetization is M ac =Xzflac- We estimate M. ic /M s m 
0.1 for a system with linear size 100 jjm for HoNi2B2C. Thus 
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the effect of ac magnetic field due to the ac current does not af- 
fect the orientations of the magnetic moments. The ferromag- 
netic resonance by is avoided because a>j < a> g . The para- 



metric resonanceS3G2l below u> g can be avoided with a small 
ff ac « M s . 

In the presence of pinning centers, it requires a critical cur- 
rent density to depin the vortex for a dc current. The ac force 
due to the ac current helps vortices to depin from the pinning 
centers, similar to thermal noiseP2H2l For Pb-In and Nb-Ta 
alloys^, it was found experimentally that the pinning effect 
becomes very weak for an ac current with frequency above 
a> c = 10 MHz. Thus the pinning effect can be greatly sup- 
pressed with the ac current with frequencies a>j » a> c used in 
our proposal. 

When the motion of vortex lattice experiences a quenched 
pinning center, it can be effectively modelled as a particle 
moving in a periodic potential. When an ac current is ap- 
plied in additional to the dc current, there will be current steps 
known as the Shapiro steps in the IV curv d 23 l 24 l The Shapiro 
steps occur when the frequency of the ac current matches 
the washboard potential for vortex due to the pinning center, 
nu>i = G • v,; c . The Shapiro steps thus can be easily distin- 
guished from the resonance with the spin-wave excitations, 



see Eq. ( 14 1. 

The current in Eq. (JTTJ does not depend on the phase differ- 
ence between the ac current and spin wave </>, thus precludes 



the existence of the Shapiro steps due to the magnetic mo- 
ments. The reasons for the absence of the Shapiro steps are as 
follows. First, the lattice constant of spin subsystem is much 
smaller than Ai, thus the vortex lattice does not feel the wash- 
board potential due to the magnetic moments. Secondly, the 
response of magnetic moments to the magnetic field associ- 
ated with the vortex lattice is linear, thus no dc current is in- 
duced in this linear region. 

In summary, we have proposed to measure the spectrum 
of magnetic excitation in magnetic materials using motion of 
vortex lattice driven by both ac and dc current in superconduc- 
tors. The resonance between the motion of vortex lattice and 
spin wave manifests as an enhancement of current at a given 
voltage. Thus the frequency of the magnetic excitation with 
the wave number determined by vortex lattice constant can be 
extracted from transport measurement. By changing the lat- 
tice constant through applied magnetic fields, the spectrum of 
the magnetic excitation up to a wave vector of order 10 nirT 1 
can be obtained. 
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